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This  report  is  a brief  survey  of  rocket  nozzle  damping  which  can  provide  a major 
source  of  damping  in  solid  propellant  motors.  In  a stability  analysis  it  would  be 
very  advantageous  to  be  able  to  assess  nozzle  damping.  Unfortunately  the  state  of 
the  art  precludes  a reliable  evaluation  for  many  types  of  solid  propellant  motors. 

Nozzle  damping  theory  was  originally  developed  for  liquid  rockets.  Typically,  there 
is  an  injector  at  the  head  end  of  a cylindrical  combustion  chamber  terminated  by  a 
nozzle,  a geometry  far  simpler  than  that  for  a typical  solid  propellant  rocket.  The 
damping  associated  with  a rocket  nozzle,  in  many  cases,  reduces  to  evaluation  of  the 
acoustic  admittance  at  the  nozzle  entrance  plane.  For  small  perturbations  the  sonic 
point  in  the  nozzle  is  unique  in  that  downstream  of  that  point  no  wave  energy  will  be 
reflected  back  to  the  chamber  since  the  wave  propagation  velocity  is  less  than  or 
equal  to  the  mean  flow  velocity.  The  analysis  thus  requires  solving  the  flow  field 
from  the  sonic  point  back  to  the  nozzle  entrance  plane.  Simpler  treatments  are 
available  for  axial  modes  when  the  convergence  section  of  the  nozzle  is  much  shorter 
than  the  wave  length  and  will  be  treated  in  a separate  section. 
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FOREWORD 


The  occurrence  of  acoustic  oscillations  in  rocket  motors  is  a 
very  troublesome  problem.  The  ultimate  goal  in  motor  stability  theory 
is  the  ability  to  predict  oscillations  in  the  design  stage  and  to  know 
what  measures  may  be  taken  to  eliminate  them  or  reduce  their  magnitude 
to  an  acceptable  level. 

Before  any  such  overall  stability  prediction  may  be  made  it  is 
first  necessary  to  understand  all  of  the  important  sources  of  acoustic 
driving  and  damping  in  a rocket  motor.  This  report  will  deal  with  one 
aspect  of  the  subject,  namely  nozzle  damping.  It  is  basically  a summary 
of  techniques  commonly  used  to  evaluate  nozzle  damping  both  theoretical 
and  experimental.  A reference  list  to  much  of  the  appropriate  litera- 
ture is  provided. 

The  work  was  performed  by  the  Naval  Weapons  Center  during  the  period 
1 April  to  20  June  1976,  was  sponsored  by  Naval  Sea  Systems  Command 
under  SEATASK  UF31-332-308 , and  coordinated  by  Naval  Ordnance  Station, 
Indian  Head,  Maryland. 

This  report  has  been  prepared  primarily  for  timely  presentation 
of  information  and  is  released  at  the  working  level. 


R.L.  DERR 

Head,  Aerothermochemistry  Division 
Research  Department 
10  August  1976 


NWC  TM  2867,  published  by  Code  608,  25  copies. 
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INTRODUCTION 


This  report  is  a brief  survey  of  rocket  nozzle  damping  which  can 
provide  a major  source  of  damping  in  solid  propellant  motors.  In  a 
stability  analysis  it  would  be  very  advantageous  to  be  able  to  assess 
nozzle  damping.  Unfortunately  the  state  of  the  art  precludes  a reliable 
evaluation  for  many  types  of  solid  propellant  motors. 

Nozzle  damping  theory  was  originally  developed  for  liquid  rockets. 
Typically,  there  is  an  injector  at  the  head  end  of  a cylindrical  combus- 
tion chamber  terminated  by  a nozzle,  a geometry  far  simpler  than  that 
for  a typical  solid  propellant  rocket.  The  damping  associated  with  a 
rocket  nozzle,  in  many  cases,  reduces  to  evaluation  of  the  acoustic 
admittance  at  the  nozzle  entrance  plane.  For  small  perturbations  the 
sonic  point  in  the  nozzle  is  unique  in  that  downstream  of  that  point  no 
wave  energy  will  be  reflected  back  to  the  chamber  since  the  wave  propa- 
gation velocity  is  less  than  or  equal  to  the  swan  flow  velocity.  The 
analysis  thus  requires  solving  the  flow  field  from  the  sonic  point  back 
to  the  nozzle  entrance  plane.  Simpler  treatments  are  available  for 
axial  modes  when  the  convergence  section  of  the  nozzle  is  much  shorter 

than  the  wave  length  and  will  be  treated  in  a separate  section. 

% 


RESUME 'OF  FUNDAMENTAL  MODELS 


A model  for  three-dimensional  waves  in  a nozzle  with  a one- 
dimensional  mean  flow  is  derived  in  Ref.  1.  It  should  be  noted  that  the 
model  is  for  a nozzle  with  a cylindrical  cross  section.*  The  one- 
dimensional  restriction  on  Man  flow  implies  that  the  nozzle  entrance 
section  is  slowly  converging.  Solutions  for  conical  nozzles  joining  a 
cylindrical  tube  are  present  - in  the  form  of  admittance  data  at  the 
nozzle  entrance.  The  authors  point  out  that  one  integration  will  result 
in  many  solutions  sines  each  integration  step  represents  a different 
motor  nozzle  geoMtry. 


*The  case  of  a two-dimensional  nossle  is  also  treated  but  will  not 
be  considered  here  and  the  discussion  will  be  limited  to  small  perturba- 
tions. 
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The  sane  basic  analysis  was  presented  in  Ref.  2.  It  differs  mainly 
in  the  fact  that  the  authors  were  able  to  show  that  an  assumption  of 
irrotational  flow  had  only  a small  effect  on  the  calculated  results.  It 
also  expanded  the  analysis  to  Include  the  solution  for  growing  (or 
decaying)  oscillations  as  well  as  for  steady-state  oscillations.  The 
report  includes  a computer  program  for  solving  the  flow  field  for  the 
convergent  section  of  a supersonic  nozzle.  The  result  of  primary  interest 
is  the  complex  admittance,  the  real  part  of  which  is  used  to  calculate 
nozzle  damping. 


EXPERIMENTAL  DETERMINATION  OF  NOZZLE  ADMITTANCE 


Crocco,  Monti,  and  Grey  (Ref.  3)  determined  nozzle  admittance  for 
axial  modes  by  direct  measurement  of  the  perturbation  velocity  and 
pressure  at  the  nozzle  plane.  The  half  power  band  width  technique  for 
axial  modes  was  used  in  Ref.  A to  measure  nozzle  damping.  Reference  5 
is  useful  in  analyzing  the  resulting  measurements  and  relating  the 
damping  results  to  the  admittance.  A third  technique  is  the  impedance 
tube  method.  This  was  the  technique  used  in  Ref.  6 to  experimentally 
determine  three-dimensional  nozzle  admittances  for  a motor  geometry  J 

typical  of  liquid  rockets.  The  experimental  results  were  found  to  be  in 
good  agreement  with  the  theoretical  results  of  Ref.  2. 

The  resultant  expressions  for  the  real  and  Imaginary  part  of 
admittance  taken  from  Ref.  6 are 

S[S2  - S2  (1  - M2)]1/2  tenh(ze)  sec2(»B)  , 

~2 = S^M 

• teah  (wa)  + tan  (*B) 

2 2 
8 + 

and 

•I*2  " <1  “ K2,,1'2  tan(»B)  sech2(«a) 

y m — ■ ^ ■ - — — 

1 IS2  ♦ (*mM)2Htenh2<»a)  + ten2(tB)J 


I 

(2) 


nozzle 


(1) 


for  axial 


modes  (8^  - 0) 


f ■ taph.Ctq)  we 

r tenh2(va)  ♦ tea2(wf) 


• IffliSll  <zq) 

teah2(«e)  4 tan2(v|) 


Equations  (3)  and  (4)  are  Identical  to  those  given  in  Ref.  7 which 
is  an  analysis  restricted  to  axial  instabilities.  In  Ref.  7 the  deriva- 
tion of  the  equation  for  the  complex  nozzle  admittance  from  which  Eq.  (3) 
and  (4)  may  be  derived  implies  that  the  wave  length  (X)  may  be  expressed 
as  2nLc>  where  n is  an  integer.  The  equation  for  X Is  given  as 

X - 2*(1  - M2) / (w/c)  (5) 

which,  for  M2  « 1,  would  be  true  for  a normal  axial  mode, 
c 


EVALUATION  OF  NOZZLE  DAMPING' 


The  equation  for  nozzle  damping  as  given  by  Ref.  8 may  be  written 
in  the  form  (Ref.  7) 

2«  fee  - - / dSjP  |2[(1  + M2)  Re{ynJ  + M + H|y2|  ] (6) 

nozzle  n 
entrance 

where  2 

f - p.  u • u, 

V - </y  dV{l/2pUl  • + 1/2  Zlj  + — Ijj  Pi  }> 

fiC  c 


and  for  the  case  of  uniform  conditions  at  the  nozzle  entrance 


2p  cVa 

n 


-[<1  + K2)  Re{yn)  + M + 5|yJ2] 


(7) 


The  above  aquations  are  baaed  on  the  analysis  of  Ref.  8. 

Using  the  aquations  derived  in  Ref.  7 the  resulting  axial  nozzle 
damping  aquation  nay  be  written  as 

* 

!sic  _ _ H+  K1  ♦ M2)/21  tanh(2wa)  (g) 

c 1 + M tanh(2ta) 

equation  (8)  is  not  applicable  for  a sabnerged  nozzle  and  is  moat 
applicable  for  the  use  when  the  parameter  is  equal  to  ooe. 
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SHORT  NOZZLE  THEORY 


The  method  in  Ref.  2 and  7 may  be  uaed  for  calculating  nozzle 
admittances  for  pure  axial  modes.  However,  much  simpler  treatments  are 
available  for  "short"  nozzles  for  which  the  nozzle  length  is  much  less 
than  the  wave  length.  From  Ref.  9 the  equation  for  nozzle  damping  may 
be  expressed  as 


- Be  i-p  (i  + if*  g») 


(9) 


which,  for  y ■ 


1.4  and  M1 2 


« 1 yields  the  result 


« - - 1.2  f-  (10) 

n Lc 

Experimental  data  (Ref.  10)  gave  the  empirical  result 

a - -1.0  f-  J (11) 

n Lc 

for  resonant  and  decay  test  results  and 

a - - 0.3  f-  J (12) 

n L 


for  impedance  tube  technique.  The  latter  equation  was  for  preliminary 
results  and  the  discrepancy  was  not  explained . 

1 To  compare  the  calculated  equation  with  the  experimental,  Zlnn  ■ 

(Ref.  9)  expressed  the  Mach  number  in  terms  of  J and  substituted  the 
result  into  Eq.  (9)  yielding  the  approximate  equation 


I 


! 


v - 3 


(13) 


I 


and  for  y » 1.4,  Eq.  (12)  reduces  to 
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For  a aore  rigorous  treatment  of  axial  mode  nozzle  damping  the 
reader  is  referred  to  Ref.  7 and  11.  The  technique  employed  involves 
the  evaluation  of  the  admittance  by  an  impedance  tube  technique.  The 
results  indicate  that  the  short  nozzle  theory  tends  to  give  somewhat  low 
results  (Ref.  11). 

The  concept  of  capture  area  is  discussed  in  Ref.  5 and  9.  The 
capture  area  is,  in  general,  greater  than  the  entrance  area  of  a nozzle 
with  a cross  section  less  than  the  port  area.  In  both  references  it  is 
shown  that  .the  measured  nozzle  admittance  will  include  damping  effects 
attributable  to  both  the  nozzle  and  the  nonuniform  flow  region  upstream 
of  the  nozzle. 


SCALING 


Efforts  have  been  made  to  develop  scaling  laws  from  which  nozzle 
admittance  measurements  on  one  nozzle  may  be  correlated  with  the  admit- 
tance of  other  nozzles.  In  brief,  for  longitudinal  modes  the  admittance 
depends  on  a dimensionless  frequency  (u>rc/t),  the  axial  distribution  of 
the  steady  state  Mach  number  in  the  convergent  section  of  the  nozzle, 
and  V (Ref.  12).  In  what'  follows,  when  the  performance  of  two  nozzles 
are  compared  the  values  of  y for  each  gas  will  be  assumed  equal.  Then 
yn  for  two  nozzles  with  equal  entrance  Mach  .numbers  will  be  the  same  if 
the  nozzle  geometries  only  differ  by  a linear  "contraction"  or  "stretch- 
ing" and  if  the  dimensionless  frequencies  are  related  by  the  expression: 

■ Fj/o  in  which  o is  the  ratio  (Zj/Zj)  of  a characteristic  length. 

A second  type  of  scaling  is  that  between  cold  flow  admittance  data 
and  actual  motor  data  (Ref.  12). 

In  order  to  compare  the  admittance  values  it  is  first  necessary 
that  the  nozzles  be  geometrically  similar.  In  addition,  the  wave 
lengths  mist  be  equal.  Then  if  M2  « 1 (Ref.  12)  the  nozzle  admittance 
(T„)  of  the  actual  motor  is  obtained  from  the  cold  flow  prototype  by  the 
relation 


<Vh  ■ <V.  <15> 

v* 


In  Ref.  12  it  is  stated  that  experimental  testing  is  required  to  estab- 
lish the  validity  of  the  previous  correlation. 
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SUBMERGED  NOZZLES  AND  MULTI  NOZZLE  CONFIGURATIONS 


In  many  solid  propellent  rocket  motors  the  nozzle  is  submerged  as 
shown  schematically  In  Figure  1.  Experimental  evidence  has  been  obtained 
that  indicates  the  submerged  nozzle  provides  less  damping  then  does  a 
"conventional"  nozzle  (Ref.  13). 


MOZZLI  umUMCC  PLAN! 


FIGURE  1.  Schematic  Representation  of  a Submerged  Rocket  Nozzle. 
Chamber  cross-sectional  area  ■ Sq;  cavity  entrance  area  » Sg; 
nozzle  entrance  area  - Sg. 


Equation  (6)  was  used  for  calculating  axial  node  damping  for  the 
geoaetry  shown  in  Figure  1 with  the  exception,  that  Janardan  and  Zinn 
(Ref.  13)  feel  that  the- integration  should  be  across  the  entire  plane  of 
the  nozzle  entrance  as  indicated.  Assuming  a one-dimensional  axial  swan 
flow  and  wave  structure  the  dimensionless  nozzle  alpha  is  then  given  by 
lq.  (16) 


oJv  <5.0  - 5.0.)  + 1(1  + 5?o  + mV)/2J  tanh(2«a) 

-A-£  ",»■  • • JS-S L* (16) 

e [1  + M tanh(2«o) ) 

6 


• As  pointed  out  in  Ref.  13  when  there  is  no  cavity  flow  the  equation 
reduces  to 


T 


I 
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o L Mo  + 1(1  + M2o  )/2J  tanh(2iio) 
n c n n n n 

c [1  + M taui(2na)] 

c 


(17) 


Thus  Che  presence  of  reverse  flow_adds  two  terms  of  opposite  sign.  The 
damping  is  decreased  by  the  term  MgOs  which  is  larger  than  the  addi- 
tional damping  term  (-1/2  as  tanh (2na) . Thus  a decrease  of  damping 

is  caused  by  recessing  the  nozzle.  It  should  be  noted  that  Eq.  (16)  and 
(17)  reduce  to  Eq.  (8)  if  the  area  ratio  is  set  equal  to  one,  a result 
which  should  be  expected.  Equation  (16)  may  be  evaluated  for  a known 
geometry  and  fluid  properties  provided  a is  determined.  In  Ref.  13,  for 
the  geometry  shown  schematically  in  Figure  1 the  authors  found  that  the 
parameter  a was  nearly  independent  of  cavity  depth  and  cavity  flow  rate. 
The  value  of  the  parameter  a may  be  determined  by  the  modified  impedance 
tube  technique,  for  example.  Ref.  6 and  7.  In  Ref.  7 it  is  pointed  out 
that  nozzle  damping  can  be  effected  by  the  number  of  ports.  They  found 
that  quadruple  port  nozzles  provided  less  damping  for  axial  instabil- 
ities then  for  single  and  dual  ports  whose  "’damping  was  approximately  the 
same. 


CONCLUSIONS 


The  current  state  of  the  art  precludes  the  direct  calculation  of 
nozzle  losses  for  many  common  solid  propellant  motor-nozzle  configura- 
tions. The  calculation  for  a cylindrical  motor  terminating  with  a 
"conventional"  nozzle  (Ac/Ag  ■ 1)  with  slowly  tapering  entrance  section 
may  be  treated  reasonably  well  for  one-  or  three-dimensional  acoustic 
oscillations  in  the  nozzle.  This  configuration  would  apply  to  many 
liquid  rocket  motors  but  not,  in  general,  to  solid  propellant  motors. 
Equations  and  experimental  data  are  available  for  evaluation  of  damping 
of  axial  modes  by  short  nozzles. 

The  current  available  approach  for  evaluating  nozzle  damping  for 
motor  geometries  and  wave  structure  not  detailed  above  is  to  evaluate 
the  nozzle  plane  admittance  by  experimental  methods,  e.g. , impedance 
tube  methods,  and  use  that  result  to  calculate  nozzle  damping. 

The  problem  of  evaluating  nozzle  damping  in  solid  propellant  rockets 
Is  a difficult  and  not  fully  solved  problem.  However,  nozzle  damping 
may,  in  many  cases,  be  an  Important  factor  when  it  comes  to  using  stabil- 
ity models  as  a means  of  predicting  stability  trends  for  motors.  In 
particular,  methods  (analytical  and/or  experimental)  are  needed  for 
determining  nozzle  damping  of  transversa  and  three-dimensional  modes  of 
oscillation  in  solid  propellant  rocket  motors  having  complex  port  and 
nozzle  geometries. 
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NOMENCLATURE 


c Sonic  velocity  _ 

F Dimensionless  frequency  (urc/c) 

J Area  ratio  of  nozzle  throat  to  nozzle  entrance  plane 

Lc  Chamber  length 

M Mach  number 

P Pressure  amplitude 

Pl  Pressure  perturbation 

rc  Chamber  radius 

S Surface  area  dJ  (x) 

Sgjj|  nth  root  of  the  equation  — jj 

Sq  Nozzle  entrance  area  x 

IT  Mean  velocity 

ui  Velocity  perturbation 

V Volume 

V See  Eq.  (6) 

Y Nozzle  admittance  u^/p^ 

y Non-dimensional  admittance  pc  (u^/p^) 

a Ratio  of  amplitude  of  reflected  pressure  wave  to  amplitude  of  incident 

wave  at  the  nozzle  entrance  defined  by 
Oq  Temporal  nozzle  decay  coefficient 

8 Phase  change  between  Incident  and  reflected  pressure  waves  at  nozzle 
entrance  defined  by  ir(l  + 28) 

Y Ratio  of  specific  heats 

X Have  length  v 

p Density 

a Ratio  of  characteristic  length  of  two  nozzles 
®n  Sn^c 

®8  8g/Sg 

u Angular  frequency 

Subscripts 

c Chamber  cross-section 

1 Imaginary 

n Nozzle  entrance  cross-section 

H Value  for  "actual"  rocket  motor 

E Value  for  prototype  motor 

r Real 

s Area  at  nozzle  entrance  plane  inclusive  of  nozzle  entrance  section 

( ) Average  value 

< > Time  everage  value  over  a cycle 
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